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Abstract. We compnte the stochastic gravitational wave prodnction from Affleck-Dine 
condensate fragmentation in the early nniverse, focnsing on an effective potential with a 
logarithmic mass correction that typically arises in gravity mediated snpersymmetry break¬ 
ing scenarios. We find that a significant gravitational wave backgronnd can be generated 
when Q-balls are being formed ont of the condensate fragmentation. This gravitational wave 
backgronnd has a distinct mnlti-peak power spectrnm where the trongh is closely linked to 
the snpersymmetry breaking scale and whose freqnencies are peaked aronnd kHz for TeV 
snpersymmetry breaking. 
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1 Introduction 

In the early universe, some scalar fields may develop large vevs, forming Affleck-Dine (AD) 
condensates [1]. This is generically expected in supersymmetric (SUSY) extensions of the 
Standard Model, as there are typically plenty of flat directions where the scalar potential 
vanishes. The Minimal Supersymmetric Standard Model already has hundreds of flat direc¬ 
tions [2, 3], which leads to interesting cosmological consequences (see [4-6] for recent reviews). 

For low energy SUSY, the flat directions are lifted by soft breaking terms not far away 
from the weak scale. These soft terms are negligible in the very early universe when the 
Hubble parameter H is much larger than the weak scale. Indeed, the finite energy density of 
the universe can induce an 0{H) effective mass for flat directions [2, 7]. The induced mass 
term can sometimes be negative, which efficiently drives the held away from the zero value 
during inhation, and, once non-renormalizable terms are included, the held can be stabilized 
at a hnite value, forming the AD condensate [2, 7]. The AD condensate almost adiabatically 
follows the decrease of 77, and, when H is surpassed by the soft mass scale, the condensate 
starts oscillating around the potential’s minimum and fragments. Therefore, the scale of AD 
condensate fragmentation is closely linked to the scale of SUSY breaking. On the other hand, 
if an AD held also carries the baryon (i?)/lepton (L) number, the observed matter-antimatter 
asymmetry can be created in this process [1]. 

Q-balls are a type of non-topological soliton, which has a localized, time-dependent and 
non-dissipative held conhguration [8, 9]. They exist in a scalar held theory endowed with a 
continuous global symmetry and an “open-up” potential [8]. Due to the global symmetry, a 
Noether charge Q can be dehned, and thus a Q-ball is a concentrated lump of Q charges. 
Physically, they can form because the spherically symmetric Q-ball conhguration is energet¬ 
ically favored in such a held theory [8]. (See [10] for the existence of long-lived composite, 
charge-swapping Q-balls in the same theory where Q-balls exist.) Q-balls are relevant in the 
early universe and can be an important consequence of SUSY hat directions ([11-31] and 
see [4-6] for recent reviews). Indeed, resonance often takes place in the decay process of 
an AD condensate: Fluctuations in the AD condensate are usually rapidly amplihed by the 
condensate’s coherent oscillation, which fragments the AD condensate into Q-balls. 

Unlike photos, gravitons decouple from the hot, dense Big Bang plasma at a very high 
temperature, thus relic gravitational waves (GWs) can be a probe into the very early universe. 
Signihcant stochastic GW backgrounds can be generated in early universe phase transitions 
where resonance takes place, such as in preheating after inhation (see [32-43] and references 
therein). This also happens in AD condensate fragmentation when the condensate has become 
dominant in the universe’s energy budget before it decays [25, 26]. In some cases, spherically 
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symmetric and energetic lumps, such as oscillons and Q-balls, can emerge as the resonance 
product in the phase transition. As spherical symmetric configurations do not emit GWs, 
we can expect that, in these cases, the GW frequency bands associated with the lump’s 
characteristic sizes are relatively suppressed, so the GW power spectrum generated has a 
distinct multi-peak structure. For example, this has been shown to be the case in oscillon 
preheating after axion monodromy inflation [43]. 

In this paper, we calculate the stochastic GW background generated when an AD con¬ 
densate fragments into Q-balls. Our numerical results will be presented for the running mass 
model (see Eq. (2.1)), which typically arises in gravity mediated supersymmetry breaking 
scenarios. The GW production of this model has been investigated with a lattice compu¬ 
tation previously in Ref [26], but we re-visit it using an independent code and with much 
higher numerical resolution. We show that the GW power spectrum generated indeed has 
a multi-peak structure, with the frequency bands of the troughs linked to the characteristic 
sizes of the Q-balls, analogous to that of oscillon preheating in axion monodromy inflation 
[43]. The peak frequency of the GW power spectrum is also more carefully determined. 


2 Model and implementation 


There are typically many flat directions in supersymmetric extensions of the Standard Model, 
but the couplings between the flat direction fields themselves as well as the couplings of the 
flat direction fields to other light fields are usually suppressed [1, 15]. Therefore, one may 
simply focus on one flat direction, which we will label as a complex scalar field (with a 
canonical kinetic term). In gravity mediated SUSY breaking, a typical effective potential for 
is given by [6, 15] 


V = + Kml/ 2 \^f log 


|( J )|2 


+ Am 


3/2 


-|- c.c. I -|- A 


| 4 > 

M; 


2n-2 


2n-6 ’ 


( 2 . 1 ) 


where K,A,X are dimensionless constants, n is an integer, Mp is the reduced Planck mass 
and c.c. stands for complex conjugate. In this scenario, the flat direction (or sfermion) mass 
is linked to the mass scale of the gravitino. Gonsidering the bounds from collider physics 
and naturalness, this mass may not be much far away from OiAoV). The logarithmic term 
parameterizes the running of the flat direction mass and mainly comes from the contribution 
of one loop gaugino corrections, and typically the value of K is roughly where 

mi /2 and mj are the gaugino and slepton masses respectively and ag is the SU{3)c coupling 
constant [4-6]. The A term and A term come from non-renormalizable interactions, with A 
and A usually assumed to be order one. In the MSSM, n is bigger than 3 and can be as big as 
9 [3], or even bigger if further symmetries are assumed [2]. As mentioned in the introduction, 
the finite energy density in the early universe induces extra Hubble parameter dependent 
terms in the potential: chH^\^\^ and /Mp~^ + c.c.'), cp and ap being constants 

[2, 7]. The cp and A terms set the initial homogeneous value for ~ Mp(m3/2/Mp)3/(’^“2)^ 
In the context of SUSY models, the U{1) asymmetry generated are typically very small 
compared to the condensate itself [2, 7] but still typically bigger than the required baryon 
number density today np/n.y ~ 10“®. This can easily be solved by one of the following 
mechanisms [25, 26]: The U{1) asymmetric terms are somehow suppressed; Or, there are 


- 2 - 





substantial entropy productions at later stages; Or, the flat direction is that of B — L = 0 
and the AD condensate decays above the scale where the active sphaleron process can destroy 
possible big B + L] Or, the flat direction is that of -B = 0 and L = 0 [19, 22, 26]. GWs are 
only significantly produced when the AD condensate starts to fragment, after which time the 
A, A, ch and an terms become negligibly small. On the other hand, the GW production only 
depends on the energy momentum tensor of the condensate and is “blind” to the matter-anti- 
matter asymmetry in the AD condensate. So for our purposes we may safely neglect these 
terms, which we have checked numerically. It is in this approximation that a global U{1) 
symmetry emerges and Q-balls can be exactly defined. 

Intuitively, Q-balls can form in this model because the AD condensate develops a neg¬ 
ative pressure when the field oscillates around the minimum [26]. In other words, there is a 
process of parametric resonance, which quasi-exponentially amplifies certain wavelengths of 
fluctuations. It follows that the AD condensate can efficiently fragment into lumps, which 
then evolve to spherically symmetric Q-balls. Technically, one may check that (cf. 

Eq. (2.1) with AT < 0, A = 0 and A > 0) has a minimum at [‘hj ^ 0, thus this potential is of 
the “open-up” type that supports Q-balls [8]. 

The “open-up” condition for Q-balls to form is very similar to that of oscillons [44]. 
Oscillons [45] are also localized, spherically symmetric, time dependent field configurations, 
and can be produced from fragmentation of an inflaton condensate in preheating [46]. Un¬ 
like Q-balls, there is no Noether symmetry to guarantee their stability, but oscillons have a 
prolonged dissipation time. We will see shortly that Q-balls and oscillons play a very similar 
role in shaping the produced GW power spectrum. 

For potential (2.1), d^U/d|<I>p has a soft, logarithmic divergence at |<h| = 0. Physically, 
small |<I>| is of course outside the validity of this effective potential, as the running mass 
correction becomes bigger than the leading term for small |<h|. Fortunately, this limitation 
of the effective potential does not significantly affect our computations. Note first that the 
equations of motion do not have any divergences. Also, the typical field value during its 
evolution is far away from the minimum. We find that numerically the effective mass only 
deviates slightly from 7713 / 2 , by about 10% in the simulations presented below. In the following, 
we will refer to the effective mass as 7773 / 2 . 

Fragmentation of the AD condensate is a highly non-perturbative process and thus 
a numerical approach is needed to understand its full dynamics. Fortunately, a classical 
approximation is applicable to the fragmentation, as the occupation numbers for the relevant 
modes involved become very high once resonance kicks in. We assume the AD condensate 
becomes the dominant component of the universe before fragmentation [26]. We have run 
lattice simulations with metric perturbations included and found that the backreaction to the 
metric is as expected negligible, so it is sufficient to solve the Klein-Gordon equation in the 
Friedmann-Robertson-Walker (FRW) background: 


d^ d 1 
^ + 3F--1$ = - 




dt 


dV 

d^* 


( 2 . 2 ) 


where a is the scale factor, H is the Hubble parameter, <I>* is the complex conjugate of 
and is the 3D flat space Laplacian. The evolution of the scale factor is determined by the 
Friedmann equation 

3Mj>H^ = {p), /9= |6|2 + |V4>|2 + U, (2.3) 

where ( ) means averaging over the space. GWs are the tensor modes of the metric 

ds^ = —dt^ -I- a^{5ij + hij)dx'^dx^, 
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where hij is transverse and traceless, dihij = ha = 0. To compute the stochastic GW 
background generated from fragmentation of the AD condensate, we solve the perturbed 
Einstein equations for the tensor modes (which are structurally similar to Eq. (2.2)), 


d 1 


h - — _ T- 

- ^2 A, 
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(2.5) 


where are the transverse and traceless part of the energy momentum tensor of the AD 
condensate. The energy density of the stochastic GW background is given by 
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where V is the space volume. Practically, one may solve an equation identical to Eq. (2.5) 
except that is replaced by Tij, the energy momentum tensor without the transverse and 
traceless projection, and, whenever an output is desired, do the projection at the level of 
Eq. (2.6) in Fourier space, which is computationally much more efficient. One of the key 
observables we are interested in is the fractional GW energy density in the current epoch (we 
will refer to this quantity as the GW power spectrum in the following), which is defined by 




Pc d In /o 


h"^ dpgw 
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where h is the observational Hubble parameter today, 0 denotes the current epoch of the 
universe, /o = k/2'KaQ is the GW frequency today and = 3Mp77g is the critical energy 
density today. 

Our lattice code is based on the public code HLATTIGE [47], where the equations of 
motion are re-cast in a different form in order to make use of more accurate, stable symplectic 
integrators. The use of a large lattice is particularly important in this study, because the 
separation of the lattice points become bigger with time, while the size of a Q-ball remains 
roughly the same. We need to make sure there are sufficient points to resolve a Q-ball at the 
end of the simulation. The results will be presented using a 256^ lattice, for which the GW 
power spectrum has sufficiently converged. The physical size of the lattice should be adjusted 
such that parametric resonance is captured in the simulations, leading Q-balls to form. This 
is achieved by choosing the lattice physical separation such that Q-balls are resolved within 
the lattice box. 

7713/2 in the potential is chosen to be ITeV for the hducial model, which applies to 
Figs. 1, 2 and 3, but we will also consider 7773/2 = lOTeV, lOOTeV in Fig. 4, to see the 
dependence of the GW peak frequency on the SUSY scales. K is typically within the range 
between —0.01 and —0.1 in MSSM [15], and we choose K = —0.1 as the hducial value here. 
We hnd that the AD condensate fragments at a slower pace for a smaller K, but the total 
GW production is insensitive to K, consistent with the previous results with a smaller lattice 
(64^) in Ref [26]. The initial value of the homogeneous mode $0 is chosen as 10^® GeV, which 
corresponds to the case where the hat direction is lifted by 77 = 9 terms in Eq. (2.1). The 
physical time scales of the fragmentation is insensitive to initial 4>o, but a smaller initial 
$0 does lead to a smaller scale factor a at the onset of the fragmentation, as one would 
expect from the Friedmann equation. As the initial net charge is usually very small (A large 


4 





net charge would be difficult to satisfy the current cosmological bounds on the matter-anti- 
matter asymmetry) and the GW production is blind to the U(l) charges, we will focus on 
a vanishing initial <ho for simplicity. We have checked that adding some initial net charge 
has little impact on our main results. There are two kinds of inhomogeneities in the AD 
condensate that can be amplihed by the parametric resonance process of the fragmentation: 
1) Quantum fluctuations which exited the horizon during inflation and have re-entered the 
horizon afterwards; These inhomogeneities are roughly scale invariant, Gaussian and of order 
|<5<1>/<1>| ~ 10“^; 2) Quantum fluctuations of the AD condensate itself, which are also Gaussian 
but much smaller than the classicalized inhomogeneities from inflation. We will thus only 
include the hrst kind as the initial inhomogeneous seeds to our simulations, and consider an 
initial power spectrum for the complex held d> 

27r^ 27r^ 

1P$/$„(/^) = A|x^ = 2.4x10-9x^. (2.9) 

(In comparison, the quantum huctuation of the AD condensate itself is roughly 7 ^$/$q(A:) ~ 
1/ ) where d>o is the homogeneous mode mentioned above. This is much 

smaller at the time of fragmentation for all the relevant modes around k ~ O(TeV).) However, 
we hnd that the total GW production is actually insensitive to the amplitude of the initial 
inhomogeneities, again consistent with the previous results with a smaller lattice (64^) in Ref 
[26]. 

3 Numerical results 

For our hducial model, linear parametric resonance takes place before t 270m3^2- After 
that, the dynamics becomes highly nonlinear and the AD condensate rapidly fragments into 
Q-balls. Q-balls have become the dominant component in the universe by t ~ 520mg^2 when 
the energy stored in regions with an energy density greater than twice the average reaches 
the plateau of about 80%; See Fig. 1. This plateau is only slowly lifted to about 90% by 
the end of the simulation t ~ 2500m[^^2 when Q-balls have properly formed. However, after 
domination, Q-balls continue to evolve towards spherical symmetry and become increasingly 
massive: After the fragmentation, the energy stored in regions with a density greater than 
100 times the average increases almost linearly to about 65% at the end of our simulation. 

We can divide the AD condensate fragmentation and its GW production into four stages 
(See Figs. 2 and 3; the time scales are for the hducial model.): 

• Linear parametric resonance (t ~ 0 —?■ During this stage, the GW power 

spectrum grows steadily at the low momentum end of the spectrum. However, there 
is no multi-peak structure in the GW power spectrum yet. The peak wavenumber 
of the linear growth of the perturbations has been estimated to be around ~ 

(1 — |iF|/4) [19], which agrees with our numerical results. 

• AD condensate fragmentation {t ~ —)• This is when the huctuations 

produced by linear parametric resonance has become sufficiently signihcant and the 
dynamics starts to become highly nonlinear. The AD condensate fragments efficiently 
during this short period. The low momentum modes of the fluctuations start to re¬ 
scatter to higher momentum modes. During this rapid process, the GW power spectrum 
grows very quickly, particularly at the higher momentum end of the spectrum. Towards 
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Figure 1. Evolution of the percentage of the integrated energy in regions with an energy density 
p > 2 (blue line) and p > 100 (red line), for the fiducial model. The average energy density is 
normalized to 1. The lattice size is 256^, the initial box size equals 0.05i?“^ and the lattice spacing 
and time step are dx = lOdt ~ 0.078m^y2- 



(a) Linear perturbation 



(c) Emerging Q-balls 



(b) Fragmentation 


1 • 

1 • 

• / 

• / 

/ 


» Wx / 

TimeStep: 1126 


(d) Properly formed Q- 
balls 


Figure 2. Evolution of the energy density via volume rendering, from the same run as Fig. 1. Note 
that the average energy density is normalized to 1, and the plotted max value of the energy density 
of (a) and (b) is 5, while for (c) and (d) it is 200. Lighting effects are used so that the minute linear 
perturbations in (a) can be seen. One “TimeStep” is 250dt = can be seen from (d), a 

Q-ball is described by 0(10^) lattice points at the end of the simulation (box size: 256^). 


the end of this stage, multiple peaks and troughs start to emerge in the GW power 
spectrum. 

• Q-balls emerging {t ~ 410m“^2 520m“^2)' After the fragmentation, Q-balls emerge. 

The multiple peaks and troughs become apparent in the GW power spectrum. Most 
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Figure 3. Evolution of the GW power spectrum. The GW power spectrum grows upwards and the 
four spectra correspond to the four stages of Fig. 2 respectively. The frequency modes in the bump of 
the blue line have grown due to parametric resonance. The black line is the final stabilized GW power 
spectrum. The vertical lines (solid and dashed) are the frequencies corresponding to m. 3 / 2 , 27713 / 2 , 
87773 / 2 , 47773 / 2 , 57773 / 2 . The odd multiples of 7773/2 align with the troughs. This is analogous to the 
case of oscillon preheating [43]. The lattice size is again 256^, but a bigger box size 0.125i/“^ is used 
to resolve the low frequencies. The lattice spacing and the time step are dx = lOdt ~ 0.197777gy^2- 


gravitational waves are produced during the first three stages, and we will refer to the 
duration of the third stages as the “fragmentation time”. 

• Q-balls become very massive {t ~ 520777~,2 After emerging, Q-balls slowly evolve 
towards spherical symmetry and steadily become more massive. The first trough slowly 
becomes increasingly apparent. The black curve in Fig. 3 is plotted around four folds 
of the fragmentation time when the GW power spectrum has become very stable. 

In Fig. 3, we see that the GW power spectrum at the frequency bands corresponding 
to the odd multiples of 7773/2 are suppressed, thus forming a multi-peak power spectrum. 
This is very much like that of oscillon preheating after axion monodromy inflation [43], in 
which case the inflaton condensate fragments into oscillons. As the physical reason for this 
to happen in AD condensate fragmentation is almost identical to that in oscillon preheating, 
we shall only mention the essential points in the following (a detailed discussion of the com¬ 
pletely analogous oscillon case can be found in [43]). Q-balls, like oscillons, formed out of a 
condensate’s fragmentation, will evolve towards spherical symmetric lumps and spherically 
symmetric lumps do not emit GWs, so the GW frequency bands associated with the char¬ 
acteristic structures of Q-balls are suppressed in the spectrum. The characteristic structures 
of Q-balls are linked to the mass parameter, so Q-balls suppress the GW production at the 
frequency bands associated with the multiples of the mass. The reason why only frequencies 
linked to the odd multiples of the mass parameter are suppressed is due to the fact that 
the potential is invariant under 4> —?■ — [43]. However, the frequencies associated with the 

higher multiples are only suppressed during early times of Q-ball formation, since the higher 
frequency spherical structures decay at late times of the formation. So we end up with a GW 
power spectrum where the frequency band associated with the mass parameter itself (i.e., the 
hrst vertical line to the left, which corresponds to the overall size of the Q-ball) is suppressed. 
For the logarithmic Q-ball case at hand, the hrst trough linked to the Q-ball size appears 
relatively late and is relatively shallow, due to the softness of the logarithmic interaction. 
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Figure 4. Final GW power spectra for 7713/2 = ITeV, lOTeV and 100 TeV. The parameters and 
initial conditions are the same as the fiducial model, except for 7713 / 2 . 


The frequencies of the GW power spectrum (as observed today) have also been carefully 
determined in our simulations, which peaks around a few kHz for the case of 7713/2 = 1 TeV 
(See Fig. 3). Thus, the stochastic GW backgrounds for this scenario can not be observed by 
the current and upcoming interferometer based GW experiments (such as advanced LIGO) 
contrary to the claims of [26]. However, new methods of observing ultra high frequency GWs 
are being developed, for example [48]. Since a bigger 777 - 3/2 leads to higher GW frequencies, 
in the event of possible observations, one may use this stochastic GW background as a probe 
to the SUSY energy scale. See Fig. 4 for the GW peak frequency’s dependence on the SUSY 
scale 7773 / 2 . 


4 Conclusion 

In summary, we have studied the dynamics of AD condensate fragmentation in the early 
universe and the associated stochastic GW production, based on the running mass potential 
(2.1) that typically arises from SUSY flat directions. Q-balls can be copiously produced in the 
AD condensate fragmentation, and a significant stochastic GW background can be produced 
when Q-balls are being formed out of the AD condensate. We have identified four distinct 
stages of this process. The generated GW power spectrum possesses a distinct multi-peak 
structure, which was unidentified in the previous studies. The trough in the GW multi¬ 
peak spectrum is linked to the SUSY breaking scale. This GW background, however, is not 
observable in the current or upcoming GW experiments (such as advanced LIGO), due to 
high peak frequencies and relatively small amplitudes. If this stochastic GW background 
can be observed in future high frequency GW detectors, one may use it as an independent 
probe to the SUSY energy scale. The SUSY energy scale this method can probe depends 
on how high the GW detector can probe a stochastic GW background’s frequency. The 
reason for the trough to form in the GW spectrum is because spherically symmetric Q-balls 
are suppressing the GW production at the relevant frequencies, analogous to the scenario of 
oscillon preheating after axion monodromy inflation. 
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